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Preface 


The  purpose  of  this  thesis  is  to  evaluate  three  possible 
choices  of  the  atomic  Hamiltonian  as  the  energy  operator  for 
an  atom  in  an  electromagnetic  field.  Confusion  between  the 
Hamiltonian  and  the  energy  operator,  and  misinterpretation 
of  operators  in  different  representations  have  plagued  stu¬ 
dent  and  author  alike.  Hopefully,  this  thesis  will  provide  a 
document  that  can  be  used  by  others  looking  at  the  atomic 
energy  operator  or  at  various  forms  of  the  Hamiltonian,  since 
there  appears  to  be  some  confusion  in  the  current  literature. 
There  are  possible  follow  on  studies  suggested  by  the  paper, 
primarily  in  experimental  confirmation  of  the  energy  operator. 
The  reader  is  expected  to  have  had  at  least  a  first  course 
in  quantum  mechanics . 

I  wish  to  sincerely  thank  Dr.  Richard  J.  Cook  for  his 
patience  and  assistance  during  the  completion  of  this  thesis, 
and  for  originally  suggesting  the  topic.  I  must  complement 
him  on  his  thorough  comprehension  of  the  subject,  a  portion 
of  which  he  has  passed  on  to  me.  Also,  I  extend  my  appre¬ 
ciation  to  Dr.  Charles  H.  Townes  for  suggesting  that  a  charged 
particle  reacts  equally  to  a  longitudinal  or  transverse  field. 
Further,  it  is  difficult  to  sufficiently  commend  my  wife 
Barbara  for  her  support  during  the  course  of  this  program, 
and  her  patience  during  the  long  hours  of  study.  Finally,  I 
wish  to  thank  Cascade  for  her  understanding  and  constant 
companionship  during  my  studies. 
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Abstract 

Various  forms  of  the  atomic  Hamiltonian  are  examined  as 
the  energy  operator  for  an  atom  in  an  electromagnetic  field, 
using  the  semiclassical  approach.  The  unperturbed  Hamiltonian 
and  the  full  atomic  Hamiltonian  have  been  favorite  choices  in 
the  past,  but  these  give  different  predictions  if  treated  as 

_  A  i  _ 

the  energy  operator.  The  A»p  versus  r*E  controversy  is 
also  examined  and  clarified.  Both  choices  have  conflicts 
with  observations  or  physical  laws.  Gauge  invariance  and  the 
Stark  effect  are  considered.  The  multipole  Hamiltonian  is 
presented  and  argued  as  the  correct  energy  operator.  Calcula¬ 
tions  for  three  photon  absorption  in  a  two  level  atom,  using 
time  dependent  perturbation  theory,  yield  significant  differ¬ 
ences  in  the  predictions  for  absorption  rate  and  resonant 
frequency. 


QUANTUM  MECHANICAL  ATOMIC  ENERGY  OPERATOR 


I .  Introduction 


Background 

In  quantum  mechanics,  the  state  vector  |  vp>  describes  a 
physical  system,  and  contains  all  information  that  can  be 
determined  about  that  system.  If  the  state  vector  is  known 
for  all  time,  then  one  can  easily  predict  or  explain  observed 
physical  phenomenon.  A  fundamental  postulate  of  quantum 
mechanics  is  that  the  Schrodinger  equation, 

ih—^l  4)>  =  H  ( t)  |i|»  (1.1) 

4  • 

governs  the  time  evolution  of  the  state  vector,  where  H(t) 
is  the  Hamiltonian  operator  for  the  system.  The  solution  of 

this  equation  provides  much  insight  for  atomic  and  molecular 
processes,  particularly  interactions  with  light.  Unfortu¬ 
nately,  the  Schrodinger  equation  is  difficult  to  solve,  and 
has  exact  solutions  for  only  a  few  simple  cases.  Generally, 
approximation  methods,  such  as  perturbation  theory,  are  used 
to  arrive  at  a  solution  which  gives  a  satisfactory  explana¬ 
tion  to  observations. 

In  the  Schrodinger  equation,  the  Hamiltonian  generally 


represents  the  total  energy  of  the  system.  The  exact  quantum 
mechanical  form  usually  originates  from  the  classical  Hamil¬ 
tonian.  In  atomic  systems,  the  form  of  the  Hamiltonian  is 


important  for  predictions,  often  through  perturbation  theory, 
of  transition  rates,  lineshapes,  multiphoton  processes,  and 
Stark  shifts.  Different  choices  for  the  Hamiltonian  can 
result  in  significantly  different  predictions.  The  key 
question  is  which  Hamiltonian  actually  represents  the  energy 
of  the  atom.  This  has  been  the  source  of  a  great  debate  and 
misconception  within  the  scientific  community. 

Some  confusion  arises  by  the  common  treatment  of  the 
Hamiltonian  and  the  energy  operator  as  equivalent.  In  fact, 
these  quantities  are  very  different  in  quantum  mechanics 
(1:326;  2).  The  energy  operator  represents  an  observable,  so 
its  expectation  value  is  the  same  in  any  representation.  In 
general,  the  Hamiltonian  is  not  an  observable.  Under  a 
transformation  to  another  representation,  it  can  describe 
something  other  than  the  energy.  This  subtle  difference  is 
important  since  the  eigenstates  of  the  energy  operator 
represent  the  stationary  states  of  the  atom.  Care  must  be 
taken  when  defining  what  form  of  the  Hamiltonian  is  to  be  the 
energy,  a  physical  quantity,  in  a  particular  representation. 
Its  expectation  value  must  be  invariant  under  transformation. 

For  the  interaction  of  light  with  matter,  a  central 
issue  is  which  form  of  the  Hamiltonian  represents  the  true 
energy  of  the  atom  (or  what  is  the  form  of  the  atomic  energy 
operator) .  Historically,  there  have  been  two  choices  for 
the  energy  operator  in  the  semiclassical  treatment.  Semi- 
classical  refers  to  quantizing  the  atom,  but  treating  the 
electromagnetic  field  classically.  The  unperturbed  atomic 

2 


vrfv  . 


I 

1 


Hamiltonian, 
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is  most  often  used  to  represent  the  atom's  energy,  where 

A 

p  is  the  canonical  momentum  operator  and  V  is  the  coulomb 
potential  due  to  the  nucleus.  Others  use  the  full  atomic 
Hamiltonian, 
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+  V  (r) 
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where  A  is  the  vector  potential  and  V  is  the  total  scalar 
potential.  Each  of  these  forms  results  in  a  different  inter¬ 
action  Hamiltonian  when  the  atom  is  in  an  electromagnetic 
field.  The  first  results  in  the  minimal  coupling  Hamiltonian, 
which,  in  perturbation  theory,  has  transition  matrix  elements 

A 

of  the  form  <n|-^A*p|m>  ,  between  states  of  the  unper¬ 

turbed  atom.  The  latter  yields  the  electric  dipole  Hamil- 

A 

tonian,  which  has  matrix  elements  <n|-er*E|m>  .  The 

difference  has  caused  a  great  deal  of  confusion  as  to  which 
one  to  use  or  if  there  is  indeed  a  difference. 

For  several  years,  the  two  forms  were  considered  equiva¬ 
lent  since  they  are  related  to  one  another  through  a  unitary 
transformation  (3).  Lamb  (4:268),  however,  noted  that  it  is 
important  to  choose  the  correct  form  since  there  is  a  differ¬ 
ence  in  the  lineshape  predictions.  Because  of  the  gauge 

A  ^  ^ 

dependence  of  A*p  ,  he  favors  the  r*E  form.  On  the 


k 
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other  hand,  textbooks  routinely  expand  the  full  Hamiltonian, 

A 

Eq  (1.3),  and  treat  the  A*p  term  as  a  pertubation  (5:459; 
6:400-401).  It  seems  well  stated  by  Worlock  (7:1327)  that, 
"the  choice  is  largely  dictated  by  taste  and  convenience." 

A 

The  r*E  form  is  found  to  be  more  useful  in  atomic  and 

A 

molecular  applications,  where  the  A*p  form  is  used  for 
energy  bands  in  crystals.  For  single  photon  resonant  ab¬ 
sorption,  there  is  essentially  no  difference  in  the  predicted 
rates;  however,  the  two  give  different  predictions  for  two- 
photon  absorption  (7;  8).  Although  it  is  evident  that 
distinct  differences  exist,  the  primary  matter  of  which  equa¬ 
tion  represents  the  energy  of  the  atom  is  still  unresolved. 

A  further  point  of  confusion  is  the  apparent  equivalence 
of  the  two  forms  under  a  unitary  transformation.  Several 

_  A.  A 

authors  have  shown  that  <n|-  — A*p|m>  and  <n|-er-E|m> 
will  yield  the  same  results  if  properly  interpreted  (9;  10; 
11) .  Under  a  unitary  transformation,  the  basis  states  of  a 
given  representation  change.  Thus,  even  though  the  eigen¬ 
states  in  each  case  have  the  same  form,  they  do  not  represent 
the  same  state.  With  this  interpretation,  the  two  forms  are 
indeed  equivalent.  If  the  states  were  considered  the  same, 
the  gauge  dependence  of  A  would  cause  the  matrix  to  be 
gauge  dependent.  A  change  of  gauge  manifests  itself  as  a 
unitary  transformation,  in  such  a  way  that  the  matrix  is 
invariant,  and  the  two  forms  are  equivalent  (12).  While 
these  arguments  are  essentially  correct,  they  do  not  fully 
answer  which  form  of  the  Hamiltonian  represents  the  energy 


of  the  atom 


Until  the  atomic  energy  operator  is  correctly  identified, 
misinterpretation  and  confusion  will  face  those  studying  or 
aPPlying  quantum  mechanics  to  light  and  matter  interactions. 
Journal  articles  still  appear  that  have  varying  interpreta¬ 
tions  of  the  energy  operator  (8;  13;  14).  Perhaps  in  pre¬ 
vious  work,  where  the  field  intensities  were  relatively  small, 
the  differences  were  unimportant.  Now,  with  high  power  lasers 
available  that  can  produce  strong  fields,  the  differences  will 
have  a  much  greater  impact. 

ij  Thesis  Objective 

There  is  an  obvious  need  for  clarification  on  the  exact 
form  of  the  energy  operator.  Accurate  predictions  for 
several  atomic  processes  depend  on  the  correct  choice.  The 
primary  objective  of  this  thesis  is  to  examine  analytically 
possible  choices  for  the  energy  operator  of  an  atom  in  a 
time  varying  electromagnetic  field.  A  second  objective  is  to 

A 

clarify  the  controversy  over  the  equivalence  of  the  A-p 
and  r*E  forms  of  the  transition  matrix  elements.  The 
second  objective  comes  as  a  natural  prerequisite  of  the  first. 
The  choices  for  the  energy  operator  include  the  two  already 
introduced,  Eqs  (1.2)  and  (1.3),  but  because  of  certain 
shortcomings  of  each,  a  third  choice  is  presented.  This  is 
the  multipole  expansion  form  of  the  Hamiltonian, 


where  the  terms  beyond  the  unperturbed  part  are  the  energy  of 
an  electric  and  magnetic  multipole.  This  is  an  infinite 
series,  but  for  practical  applications,  it  is  truncated  to  a 


few  terms.  Here  only  the  electric  dipole  term  is  considered, 
but  in  those  instances  where  the  magnetic  dipole  or  electric 
quadrapole  is  important,  the  corresponding  terms  should  be 
included.  Treatment  of  the  multipole  expansion  form  as  the 
atomic  energy  operator  has  as  of  yet  been  unexplored. 

It  is  quite  clear  that  due  to  the  extensive  nature  of 
the  controversy,  the  only  true  test  of  a  particular  form  of 
the  energy  operator  is  ultimately  to  be  confirmation  by 
experiment.  As  part  of  arguments  for  the  multipole  form, 
this  thesis  details  calculations  for  a  three-photon  absorp¬ 
tion  experiment.  There  are^signif icant  differences  between 
the  predictions  of  the  three  possible  energy  operators. 

Final  proof  of  this  thesis  will  have  to  wait  for  such  an 
experiment  to  be  performed,  but  the  objective  of  examining 
possible  choices  of  the  atomic  energy  operator  is  met. 

Assumptions 

The  thesis  considers  the  semiclassical  approach  in  the 
analysis,  for  reasons  of  simplicity.  The  full  quantum 
approach  is  deemed  unnecessary  for  the  arguments  presented, 
and  would  significantly  increase  the  complexity.  In  addition 
only  the  electric  dipole  term  of  the  multipole  expansion  is 
considered.  The  magnetic  dipole  or  electric  quadrapole  terms 
could  easily  be  incorporated  into  the  model,  if  necessary. 


For  ease  of  notation,  the  bra-ket  convention  is  used 
throughout . 


The  experimental  calculations  use  a  general  two  level 
system  as  its  model.  Several  physical  systems  can  be  approx¬ 
imated  as  a  two  level  system,  with  good  results.  The  per¬ 
turbation  calculations  become  manageable  with  this  assumption. 
Finally,  in  the  calculations,  the  effects  of  spontaneous 
emission  is  ignored.  Again,  this  serves  to  reduce  the  diffi¬ 
culty  without  hampering  the  arguments. 

General  Approach 

The  second  thesis  objective  is  to  clarify  the  contro¬ 
versy  between  the  unperturbed  Hamiltonian,  Eq  (1.2) ,  and  the 
full  atomic  Hamiltonian,  Eq  (1.3).  Chapter  II  discusses  the 
origin  and  application  of  the  two  Hamiltonians.  Emphasis  is 
on  the  different  set  of  basis  states  each  defines  when  con¬ 
sidered  as  the  energy  operator.  The  matter  of  the  equivalence 
of  the  transition  matrices  is  also  covered  here. 

Chapter  III  addresses  the  shortcomings  of  the  first  two 
choices  for  the  energy  operator.  The  true  energy  operator 
must  be  gauge  invariant  and  must  account  for  the  Stark  effect 
for  both  longitudinal  and  transverse  fields.  Considering  the 
unperturbed  Hamiltonian  as  the  energy  operator  results  in 
gauge  dependent  predictions,  and  it  does  not  show  a  Stark 
shift.  The  full  atomic  Hamiltonian  has  a  Stark  shift  only 
for  a  longitudinal  field,  which  violates  special  relativity. 


The  multipolar  Hamiltonian  is  presented  in  Chapter  IV. 
The  justification  is  based  on  the  classical  multipole  poten¬ 
tial  and  that  the  Hamiltonian  produces  the  correct  Lorentz 
force  in  the  classical  limit.  The  Hamiltonian  is  gauge  in¬ 
variant  and  has  a  Stark  shift  for  the  total  field.  Finally, 
the  transition  matrix  elements  have  a  much  different  form 
than  the  previous  cases,  which  should  result  in  new  predic¬ 
tions  . 

Chapter  V  contains  a  perturbation  analysis  of  the  three 
cases  for  three  photon  absorption.  The  calculations  show 
order  of  magnitude  differences  in  the  absolute  rates,  and 
shifted  resonant  frequencies.  These  results  provide  a 
possible  means  to  experimentally  confirm  the  correct  energy 
operator. 


•  • 
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II .  Historical  Controversy 

Origin  of  the  Hamiltonians 

Classically,  the  Hamiltonian  of  an  isolated  single  elec¬ 
tron  atom  is  represented  by  the  sum  of  the  electron's  kinetic 
and  potential  energy.  In  the  quantum  formalism,  the  physical 
quantities  are  replaced  by  operators  to  arrive  at 

K  -  is  +  v(?)  (2*D 

the  unperturbed  Hamiltonian,  p  is  the  canonical  momentum 

operator,  -  iftV  ,  and  V  is  the  coulomb  potential  due  to  the 
e2 

nucleus,  -  —  .  Hq  defines  a  complete  set  of  basis  states 

through  the  time  independent  Schrodinger  equation, 

H0ln>  =  En | n>  (2.2) 

| n>  are  the  basis  vectors  which  can  describe  the  state 
vector : 


U>  =  ^  an  !  n>  (2*3) 

n 

Each  of  the  basis  states  is  also  called  a  stationary  state, 
and  En  represents  the  energy  of  that  state.  Most  authors 
use  these  basis  states  for  calculations  of  transition  rates, 
and  E^  as  the  energy  levels  of  the  atom.  This  method,  which 
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By  using  operators  in  place  of  the  physical  quantities,  the 
quantum  mechanical  form  of  the  total  Hamiltonian  emerges: 


+  e<t>(r,t) 


(2.8) 


e 

c 


~  —  \  ^ 
A(r, t) j 


In  the  semiclassical  treatment,  the  field  variables  remain  as 
C-numbers.  The  scalar  potential  here  is  the  total  scalar 
potential,  which  includes  the  coulomb  potential  of  the 
nucleus . 

The  primary  justification  for  using  Eq  (2.8)  as  the 
total  Hamiltonian  is  that  the  equation  of  motion  for  the  posi¬ 
tion  operator  has  a  form  that  corresponds  to  the  classical 
Lorentz  force  (6:178-179).  Using  the  Heisenberg  equation  of 
motion, 


(2.9a) 


(2.9b) 


The  brackets  refer  to  commutation.  From  these  relations,  in 
a  non-trivial  derivation,  the  results  become 


where 


eE  + 


1  e  /dr 

2  c  I  dt 


x  B 


(2.10) 
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cat 


(2.11a) 
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B  =  V  X  A 


(2.11b) 


Perturbation  Theory 

In  most  applications,  the  Schrodinger  equation  is  impos¬ 
sible  to  solve  exactly,  so  approximation  methods  are  often 
used.  A  common  method  for  time  dependent  Hamiltonians  is 
perturbation  theory,  in  which  the  total  Hamiltonian  is  put  in 
the  form  of  an  unperturbed  Hamiltonian  and  a  perturbation 
term: 


H  =  Hq  +  HI  (2.12) 

The  unperturbed  Hamiltonian  satisfies  the  time  independent  ( 
Schrodinger  equation,  resulting  in  a  set  of  basis  states. 

Ho|n>  =  En|n>  (2.13) 

The  set  of  basis  states  represents  the  stationary  states  of 

A 

the  atom,  and  depends  heavily  on  the  choice  for  HQ.  The 
state  vector  or  wave  function  is  then  expanded  in  terms  of 
these  basis  states: 

4-  =  J^an(t)|n>  e“lEnt/'h  (2.14) 

n 

The  time  dependent  coefficients,  an(t)  ,  are  interpreted  such 
2 

that  P  =  | a  |  is  the  probability  of  being  in  state  | n> 


Substituting  Eq  (2.14)  into  the  Schrodinger  equation, 
(1.1),  and  making  use  of  Eq  (2.13), 

Vfc)  =  TH2Z  <n I |m>  am(t)eiu,nmt  (2.15) 

m 

where  »  (En  -  Em)  /h  .  These  equations  define  the  time 

evolution  of  the  state  vector  in  terms  of  the  a  coeffi- 

n 

cients.  The  matrix  element  <n|HI|m>  represents  the  trans¬ 
ition  from  state  |m>  to  state  jn> 

The  transition  process  depends  heavily  on  the  interaction 
Hamiltonian  and  the  precise  eigenstates  chosen  as  the  basis. 
Special  care  must  be  taken  in  this  latter  point.  Ideally, 
and  what  is  most  often  assumed,  is  that  the  basis  states 
represent  the  true  energy  eigenstates  of  the  atom,  and  the 
eigenvalues  En  are  the  true  energy  of  the  states.  The 

A 

Hamiltonian  chosen  to  represent  the  energy  of  the  atom,  Hq 
in  this  case,  determines  the  basis  states.  It  also  determines 
the  form  of  the  interaction  Hamiltonian.  The  scientific 
community  appears  to  be  split  over  which  form  of  the  Hamil¬ 
tonian  to  use  as  the  energy  operator.  The  two  major  choices 

A 

are  the  unperturbed  Hamiltonian  Hq  ,  Eq  (2.1),  and  the  full 

A 

atomic  Hamiltonian  H  ,  Eq  (2.8).  They  each  yield  a  differ¬ 
ent  set  of  eigenstates  and  a  different  interaction  Hamil¬ 


tonian. 


If  the  unperturbed  Hamiltonian  is  the  energy  operator, 


then  the  energy  eigenstates  satisfy 

A 

H  I n>  =  E  I n> 
o'  n ' 


(2.16) 


Expanding  the  total  Hamiltonian,  Eq  (2.8),  for  an  atom  in  an 
electromagnetic  field. 


A 

H 


e<t> 


(2.17) 


assuming  that  A  can  be  considered  constant  over  the  extent 

A 

of  the  atom,  so  that  A  and  p  commute.  This  is  the  electric 

_  • 

dipole  approximation.  Note  that  A(t)  still  varies  with 
time,  and  the  potential  <j>  is  the  coulomb  potential  of  the 
nucleus.  The  interaction  Hamiltonian  is  then 


A 


(2.18) 


The  expanded  form  of  the  Schrodinger  Eq  (2.15)  becomes 

i  A  2  -  •  (jj  . 

a(t)  =  -rrr  /  <n  I  -  —  A  •  p  +  e-  A^|m>a  e1  mn  (2.19) 

in  i—j  '  me  - 2  1  m 

m  2mc 

2 

The  A  term  in  the  transition  matrix  adds  a  common  phase  to 
the  amplitudes  an  ;  thus,  it  can  be  ignored.  Most  texts 

A 

disregard  the  term  since  it  is  small  compared  to  the  A  •  p 


J 

i 

* 


ii 


I  4  • 

H 


term  (6:400).  The  transition  matrix  reduces  to 

A 

<n  I  -  —  A  •  p  I  m>  ;  this  is  the  minimal  coupling  inter- 
mc  ' 

action  form. 

A 

The  A  •  p  form  of  the  transition  matrix  is  commonly 
used  in  calculations  to  predict  transition  rates,  lineshapes, 
etc.  The  basis  states,  ln>  ,  are  assumed  to  be  the  eigen- 

A 

states  of  the  energy  operator,  Hq  in  this  case,  and  the 
transitions  are  between  these  states.  However,  one  obvious 
issue  arises.  A  is  the  vector  potential,  which  is  not 
unique  for  a  given  electromagnetic  field.  A  different  choice 
of  potentials, 

a"  =  A  +  VG  *  '  “  $  “  c  3t  (2.20) 

where  G(r,t)  is  an  arbitrary  function,  yield  the  same  E 
and  B  fields,  according  to  Eq  (2.11).  These  gauge  trans¬ 
formations  do  not  alter  physical  quantities,  but  transforming 
A  ,  in  this  case,  changes  the  value  of  the  transition  matrix. 
The  interpretation  of  the  matrix  elements  as  representing 
transitions  between  the  energy  eigenstates  results  in  gauge 
dependent  predictions,  although  this  is  a  common  practice  in 
calculations . 


Si 


k 


► 

L- - 


H  as  the  Energy  Operator 

Another  choice  for  the  energy  operator  is  the  full  atomic 
Hamiltonian,  Eq  (2.8).  This  defines  a  different  set  of  basis 
states  which  are  time  dependent,  since  the  Hamiltonian  is 


15 


time  dependent: 


H  ( t)  |  n  ( t)  >  =  En|n(t)  > 


(2.21) 


It  will  be  shown  that  the  energy  eigenvalues  are  constant  and 
equal  to  the  eigenvalues  of  the  unperturbed  Hamiltonian. 

To  justify  use  of  the  full  Hamiltonian  as  the  energy 
operator  consider  an  energy  operator  equal  to  the  sum  of  the 
potential  and  kinetic  energies 


~  i  2  — 

E  =  j  mv  +  v (r) 


(2.22) 


I  1* 


^  l  _  _ 

In  an  electromagnetic  field,  v  =  —  (p  -  —  A)  and  V(r)  =  e 4>  , 

m  e 

which  upon  substitution,  is  Eq  (2.8). 

The  eigenstates  of  the  energy  operator  are  derived  by  a 

A 

unitary  transformation,  S  ,  of  the  time  independent  Schro- 
dinger  equation.  If  H  does  indeed  represent  the  energy  of 
the  atom,  then  the  eigenvalues  are  unchanged  in  the  trans¬ 
formation,  and  Eq  (2.21)  transforms  as 


A  * 

H  |n>*  =  En|n>* 


(2.23) 


where 


A  >  AAA 

H  =  S  H  S' 


(2.24) 


n>*  =  S I n (t) > 


(2.25) 
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With  the  unitary  transformation  defined  as 


S 


ier 


A  ( t ) 


(2.26) 


then 


"*=2S+e*  (2‘27) 

A  A 

H  has  the  same  form  as  the  unperturbed  Hamiltonian,  and  has 
eigenstates  |n>*  that  are  of  the  same  form  as  the  previous 
case.  However,  the  unitary  transformation  causes  a  change  of 
representation,  so  the  new  eigenstates  are  the  same  as  the 
unperturbed  states.  The  true  energy  eigenstates  are  |n(t)>  , 
which  are  related  to  the  unperturbed  eigenstates  by 

| n  ( t ) >  =  S+ | n>*  (2.28) 


Again,  the  energy  eigenvalues  are  the  same,  which  can  cause 

A  A 

some  confusion,  but  H  and  Hq  result  in  very  different 
sets  of  energy  eigenstates. 

A 

Choice  of  H  as  the  energy  operator  also  results  in 
different  physical  predictions.  Expanding  the  state  vector 
in  terms  of  the  energy  eigenstates, 

.  E  t/h 

b(t)>  =Zwcn(t)  I  n  (t)  >  e"1  (2.29) 

m 
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Substituting  Eq  (2.29)  into  the  Schrodinger  equation. 


=  “E<n(t)llt  I  *<t)>Cm(t)e 

m 


ito 


mnt 


(2.30) 


Using  Eq  (2.28)  this  becomes 

•  i  »  1.  _  ico  , 

Cn(t)  =  Thl2*<n\  ~  er-E|m>*Cme  mnt  (2.31) 

m 

The  position  operator  has  the  same  form  in  both  representa¬ 
tions,  since  it  is  unaffected  by  the  unitary  transformation. 

/v 

The  transition  matrix  can  then  be  written  as  <n|-er«E|m>  , 

where  the  eigenstates  are  those  of  the  unperturbed  Hamil¬ 
tonian.  This  is  the  electric  dipole  interaction  form. 

The  Schrodinger  equation  for  the  two  Hamiltonians, 

Eqs  (2.19)  and  (2.31),  have  the  same  form  except  for  the 
matrix  elements,  or  more  specifically,  the  interaction  Hamil¬ 
tonians.  In  both  cases,  the  matrix  elements  are  calculated 
with  the  unperturbed  eigenstates,  even  though  they  are  differ- 
ent  than  the  energy  eigenstates  for  H  .  The  minimal  coup¬ 
ling  form  is  gauge  dependent  in  this  interpretation,  due  to 

A 

its  dependence  on  the  vector  potential  A  .  On  the  other 
hand,  the  electric  dipole  form  is  a  function  of  the  electric 
field,  which  is  gauge  invariant.  If  the  matrix  elements  are 
to  describe  some  physical  process,  the  value  cannot  depend  on 
the  gauge.  This  demonstrates  a  clear  difference  between  these 
two  forms  of  the  energy  operator.  Unfortunately,  there  is 


to 
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still  much  confusion  since  the  two  forms  are  related  by  a 


unitary  transformation.  In  quantum  mechanics,  physical  pre¬ 
dictions  are  not  altered  by  unitary  transformations;  thus, 
the  two  forms  have  been  considered  equivalent.  Yet,  it  is 
obvious  that  there  is  a  difference.  This  problem  deserves 
some  investigation  due  to  the  wide  spread  misuse  of  this 
equivalence . 

A  /V 

Equivalence  of  A_.jo  and  r  «E 

Several  authors  have  addressed  the  equivalence  of  the 

A  A 

A* p  and  r*E  forms  of  the  transition  matrix  (3;  9;  11;  12), 
but  due  to  widespread  misinterpretation  and  its  relevance  to 
this  thesis,  it  deserves  mention  here.  The  main  problem  is  in 
the  interpretation  of  the  "unperturbed"  basis  states  in  each 

<3> 

transition  matrix.  It  is  assumed  that  both  refer  to  the  same 
set  of  basis  states,  but  this  is  not  the  case.  The  numerical 
values  are  the  same,  but  they  are  states  in  two  different 
representations . 

Consider  first  the  full  atomic  Hamiltonian,  Eq  (2.8), 

H  =  ^  (P  -  ■§  A)  2  +  e<D  (2.32) 

A 

In  this  equation,  p  is  the  canonical  momentum  operator,  not 
the  kinetic  momentum.  The  kinetic  momentum  is  given  by 


mv  =  (p  -  |  A) 
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(2.33) 


Expanding  Eq  (2.32)  in  the  electric  dipole  approximation, 


e^) 


e 

me 


A  •  p 


+ 


(2.34) 


The  first  two  terms  are  in  the  form  of  Hq  ,  the  unperturbed 
Hamiltonian,  but  the  true  unperturbed  Hamiltonian  differs  in 

A 

that  the  canonical  momentum  p  is  also  the  kinetic  momentum. 

A 

Expanding  on  the  basis  states  of  Hq  as  before,  the  Schro- 
dinger  equation  becomes 


a„  =  -rr-y<n  -  —  A  •  p  m> a  e 
n  lhZ— <  1  me  r 1  m 

n 


XU!  t 

nm 


(2.35) 


The  electric  dipole  form  of  the  Hamiltonian  can  be  de¬ 
rived  by  applying  the  unitary  transformation  of  Eq  (2.26)  to 
the  Schrodinger  equation.  The  Schrodinger  equation  becomes 


ih  |^U>  '  = 


(2.36) 


where 


j>  = 


S  |  '])> 


(2.37) 


A  /  A  A  A  I 

H  =  s  h  s 


(2.38) 


Substituting  Eq  (2.32)  into  (2.38), 


H  =  E—  +e<t>-er*E 
2m 


(2.39) 


where  E  is  given  by  -1/c  3A/3t  .  Again,  the  first  two 

terms  have  the  form  of  the  unperturbed  Hamiltonian.  Expanding 

A 

on  the  basis  states  of  H  , 


i  =  ±r  V. 

n  1T1 


<n  -  er*E  m>  a  e 
t  x  1  m 


iuj  t 
nm 


(2.40) 


Thus,  it  appears  that  Eqs  (2.35)  and  (2.40)  are  identical 
except  for  the  transition  matrix.  This  misinterpretation  has 
plagued  scientists,  and  lured  them  to  believe  that  it  is 
equivalent  to  use  either  form  (7:1327). 

Upon  careful  examination,  two  mistakes  emerge.  In  the 
unperturbed  Hamiltonian  of  Eq  (2.39),  p  actually  represents 
the  kinetic  momentum,  unlike  the  earlier  case.  Consider  the 
unitary  transformation  of  the  kinetic  momentum. 


mv  =  S  mv  S+  =  S(p  -  |  A)S+ 


(2.41) 


S  and  A  commute,  but 


A  A  a  Q 

S  p  =  p  S  +  ifi  ^ 
ar 


=  p  S  +  -  A  S 


(2.42) 


s 


!• 


mv  =  p 


(2.43) 


The  basis  states  of  Eq  (2.40)  are  distinct  from  the  those  of 
Eq  (2.35)  . 

The  second  mistake  involves  the  interpretation  of  the 


coefficients  a„  .  It  was  assumed  that  in  both  cases,  a 
n  '  n 


represents  the  expansion  coefficients  of  the  state  vector  on 
the  unperturbed  basis  states.  The  transformation  of  the  state 
vector  must  also  transform  the  basis  states.  This  is  a  mis¬ 
take  often  made  in  the  literature  (13:625,631). 


A  A 
'  =  S|^>  =y  a^S | n> 


(2.44) 


Recall  Eq  (2.25) 


I  n  ( t)  >  =  S  [  n> 


(2.45) 


Thus,  in  the  second  case,  aR  are  the  expansion  coefficients 


of  the  time  dependent  basis  states. 

Eqs  (2.35)  and  (2.40)  are  indeed  related  to  one  another 
and  can  be  considered  equivalent  if  properly  handled.  Physi¬ 
cal  predictions  will  then  be  invariant  in  either  case  using 
the  correct  interpretation.  Unfortunately,  this  still  does 
not  address  the  problem  of  the  form  of  the  atomic  energy 
operator.  The  two  possibilities  still  exist;  the  unperturbed 

A  A 

Hamiltonian  H_  or  the  full  atomic  Hamiltonian  H  .  A  better 


m 


■  ■  i 
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III.  Limitations  of  and  H  as,  the  Energy  Operator 

The  correct  choice  of  the  energy  operator  must  be  consis¬ 
tent  with  experimental  observations  and  physical  laws.  The 
results  of  physical  predictions  must  be  gauge  invariant,  since 
the  choice  of  gauge  is  arbitrary.  Observed  phenomenon  such  as 
the  Stark  effect  must  also  be  accountable  in  the  theory. 
Finally,  the  theory  must  be  consistent  with  the  theory  of 

A  A 

special  relativity.  Neither  Hq  nor  H  meet  these  require¬ 
ments,  and  thus  do  not  accurately  represent  the  true  energy. 

Gauge  Invariance 

Electric  and  magnetic  fields  can  be  described  in  terms 

of  the  vector  and  scalar  potentials,  Eq  (2.11)  .  However,  the 
<8> 

potentials  are  not  unique;  a  different  set  of  potentials  from 
Eq  (2.20)  yields  the  same  fields.  This  demonstrates  that  the 
potentials  are  unphysical,  and  only  the  fields  are  important 
physically.  An  additional  condition,  a  gauge  condition,  can 
be  imposed  on  the  potentials  without  affecting  the  outcome. 

An  often  used  condition  is  the  coulomb  gauge: 

V  •  A  =  0  (3.1) 

A 

With  this  condition,  the  operators  A  and  p  commute,  and 
the  electric  field  can  conveniently  be  expressed  in  terms  of 
its  longitudinal  and  transverse  components  (16:125).  The 
longitudinal  component  of  a  vector  field,  V  ,  is  defined  as 


0 


(3.2) 


V  X  VL  = 

i 

The  transverse  component  obeys 


V  .  VT  =  0  (3.3) 

The  electric  field,  in  coulomb  gauge,  is  then  E  =  EL  +  ET  , 
where 


EL  =  -  V4> 

? _ ill 

ET  '  c  Jt  (3.4) 

The  magnetic  field  is  always  transverse: 

B  =  B T  =  7  x  A  (3.5) 

Expressing  the  full  atomic  Hamiltonian  in  this  gauge  is  use¬ 
ful,  although  not  necessary.  Any  accurate  physical  model  does 
not  depend  on  the  gauge  chosen  for  its  predictions. 

Since  the  full  atomic  Hamiltonian,  Eq  (2.8),  depends  on 
the  potentials,  its  behavior  under  gauge  transformations 
should  be  understood.  The  Schrodinger  equation  is  form  invar¬ 
iant  under  a  gauge  transformation  (1:323-325),  meaning  that 
Eq  (1.1)  takes  the  same  form  in  any  gauge: 


i*  U>  ' 


=  H  (t)  |i|» 


’  -  -  *'(t))2  +  e*'(t) 


(3.6' 


Although  p  does  not  change  form  in  the  new  gauge,  it  repre- 

A 

sents  a  different  quantity.  Remember,  p  is  the  canonical 

Ay 

momentum,  which  is  not  a  true  observable.  The  change  in  p 
results  in  a  change  in  the  basis  states  in  the  new  represents- 

A 

tion  if  the  eigenvectors  of  Hq  are  used. 

A  gauge  transformation  is  equivalent  to  a  unitary  trans¬ 
formation  to  another  representation  (1:321),  given  by 


U  =  exp 


1  t 


G(  r,t) 


(3.7) 


where  G  is  an  arbitrary  function  which  determines  the 

/\ 

gauge.  For  an  operator,  T  ,  to  represent  a  physical  quan¬ 
tity,  the  following  condition  must  be  met: 


A  <*  A  A  A  Aj 

T  =  TueU  T  U 


(3.8) 


A  __  A 

where  T  "( A  ",  $  ")  is  the  operator  in  the  new  gauge,  and 

is  the  transformed  operator.  For  example,  consider  the  true 
momentum  given  by  Eq  (2.33), 


mv  =  (p  -  —  A  ) 


(p  -  f  A  -  f  VG) 


(3.9) 
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(3.10) 
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mv 


u 


A 


As  expected,  the  true  momentum  is  gauge  invariant,  v '  =  vu 
For  a  Hamiltonian,  this  is  not  generally  true,  unless  it  is 
the  true  energy  operator  of  the  atom. 

/v 

If  hq  is  considered  as  the  energy  operator,  then  its 

A 

exact  meaning  changes  based  on  the  interpretation  of  p 
in  the  given  representation.  Even  though  the  form  remains 

/s  Jl  /s 

the  same,  HQ  is  an  unphysical  quantity,  as  is  p  .  Hq  is 
gauge  dependent  in  this  manner.  Another  way  to  see  this 
result  is  to  consider  the  transition  matrix  element, 

<n|A«p|m>  .  The  basis  states  are  the  energy  eigenstates  of 

Hq  ,  but  A  is  gauge  dependent.  Thus,  physical  predictions 

A 

become  gauge  dependent  if  Hq  is  treated  as  the  energy 
operator. 

A 

Next  consider  the  full  atomic  Hamiltonian,  H  ,  as  the 
energy  operator.  Again,  the  operator  must  be  gauge  invariant 
to  represent  a  physical  quantity.  At  first  glance,  the  opera¬ 
tor  appears  to  be  gauge  dependent  because  of  the  gauge  trans¬ 
formation  of  the  scalar  potential  in  H  : 


=  k  (p  - 1  A')2 


+  e<t>  - 


e 

c  at 


e 


(3.11) 


is  the  scalar  poten- 


i  t « 


This  is  true  only  if  the  potential,  4  , 

tial.  However,  in  the  energy  operator,  the  potential  is  ac¬ 
tually  the  coulomb  potential, 


which  is  gauge  invariant  since  charge  q  and  position  r 
are  invariant  quantities.  Hence,  the  Hamiltonian  is  also 
gauge  invariant. 

The  total  potential  can  be  expressed  as  the  sum  of  the 
potential  due  to  the  nucleus,  4>  ,  and  that  due  to  external 

sources,  4ext  / 


*  =  *0  +  W  (3*13) 

The  second  term  is  then  expanded  with  a  Taylor  series  to 
arrive  at  a  multipole  expansion, 

<t>ext  =  <t>  ( 0 )  +  V<t>{0 )  •  r  +  ... 

=  4(0)  -  EL*r  +  ...  (3.14) 

Ignoring  the  constant,  <t>(0)  ,  and  the  higher  order  terms 

A  A  /S 

H  =  k-  (p  -  f  A)2  +  e<h  -  e7-ET  (3.15) 


Thus,  the  full  atomic  Hamiltonian  is  shown  to  be  gauge 


l 


*• 
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invariant  when  handled  in  this  way.  Note  also,  that  the 

A 

transition  matrix  <n|-er»E|m>  ,  which  results  from  the  choice 

A 

of  H  ,  is  also  gauge  invariant,  since  it  depends  upon  the 
electric  field,  and  the  states  |n>  are  gauge  invariant. 

Stark  Effect 

The  Stark  effect  is  a  shift  in  the  atomic  energy  levels 
for  an  atom  in  an  external  electric  field,  and  has  been  widely 
observed  and  reported.  The  true  energy  operator  for  the  atom 
must  reflect  this  energy  shift.  The  unperturbed  Hamiltonian 
contains  no  terms  which  include  external  fields;  hence,  its 
energy  eigenstates  are  independent  of  the  external  fields. 

A 

No  Stark  effect  would  be  observed  if  HQ  were  the  atomic 
energy  operator.  The  external  field  would  simply  cause  trans¬ 
itions  between  the  eigenstates.  The  lack  of  a  Stark  shift  is 
a  strong  argument  against  the  unperturbed  Hamiltonian  as 
the  energy  operator. 

The  full  Hamiltonian,  Eq  (3.15),  includes  the  term 

A 

-er.ET  ,  so  the  energy  eigenvalues  are  a  function  of  the 
L) 

external  field.  Transforming  the  energy  operator,  Eq  (3.15) 
with  the  unitary  transformation  of  Eq  (2.26), 


A  A  A  A  i 

H  '  =  S  H  S 


Ho  _er * EL 
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^  . 


(3.16) 


The  energy  eigenvalues  of  H  "  ,  which  are  equal  to  those  of 

A 

H  ,  exhibit  a  dependence  on  the  longitudinal  portion  of  the 
electric  field.  The  eigenvalues  do  not  depend  on  the  trans¬ 
verse  component.  A  transverse  field  causes  transitions 
between  the  eigenstates.  This  conclusion  is,  however, 
erroneous . 

A  charged  particle  cannot  distinguish  between  the  longi¬ 
tudinal  and  transverse  portions  of  the  electric  field.  If 
taken  individually,  the  components  are  unphysical.  Only  the 
two  combined  together  have  a  physical  consequence.  It  would 
seem  that  the  atom  should  then  exhibit  a  Stark  effect  for 
both  the  longitudinal  and  transverse  portions  of  the  electric 
field. 

The  longitudinal  electric  field  in  coulomb  gauge  is 


Et  =  -V<J> 


(3.i7: 


1  |r-r. 


which  is  the  instantaneous  coulomb  field.  This  shows  the 
unphysical  nature  of  EL  ,  since  an  electric  field  travels 
at  the  speed  of  light,  not  instantaneously.  If  the  Stark 
effect  were  dependent  on  the  longitudinal  field  only,  then 
there  could  be  instantaneous  knowledge  of  change  in  a  charge 
distribution  at  a  large  distance,  in  violation  of  special 
relativity.  Thus,  the  Stark  effect  must  depend  on  the  total 
electric  field.  In  the  total  electric  field,  the  unphysical 
portion  of  the  transverse  component  cancels  the  unphysical 


'TTT' 


part  of  the  longitudinal  component,  insuring  that  the  field 
obeys  special  relativity  (16sl25;  17).  The  total  Hamiltonian, 
which  exhibits  a  Stark  shift  for  only  the  longitudinal  field, 
is  then  not  the  proper  choice  for  the  energy  either. 

Another  Hamiltonian 

The  two  choices  normally  considered  for  the  energy  opera¬ 
tor  have  been  shown  to  have  some  deficiencies.  A  strong  argu¬ 
ment  exists  against  the  unperturbed  Hamiltonian  because  of  its 
gauge  dependence  and  lack  of  a  Stark  shift.  The  full  Hamil¬ 
tonian,  although  gauge  invariant,  shows  a  Stark  shift  for  a 
longitudinal  electric  field  only,  which  is  in  violation  of 
special  relativity.  These  arguments  suggest  that  there  must 
be  yet  another  form  of  the  Hamiltonian  which  represents  the 
energy  of  the  atom. 

A  reasonable  choice  is  the  full  multipolar  Hamiltonian, 
which  consists  of  the  unperturbed  Hamiltonian  plus  the  multi¬ 
pole  energy  expansion.  In  the  electric  dipole  approximation, 
this  is 

A 

^  p2  -  ~  -  _ 

H  =  +  V (r)  -  er  •  (EL  +  ET)  3.18) 

This  form  reacts  to  both  the  longitudinal  and  transverse 
electric  fields  equally;  thus,  there  is  no  violation  of  spe- 

A 

cial  relativity.  In  addition,  if  V(r)  is  the  coulomb 
potential  of  the  nucleus,  the  Hamiltonian  is  also  gauge 
invariant.  Thse  topics  are  discussed  in  more  detail  in 
the  next  section. 
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I 


IV.  Multipole  Hamiltonian  as  the 


Operator 


The  Multipole  Hamiltonian 

The  multipole  form  of  the  atomic  Hamiltonian  is  usually 
derived  from  a  unitary  transformation  of  the  full  atomic 
Hamiltonian,  Eq  (2.8)  (16:166-167;  18).  The  final  result  is 


i_  r-  r1 
2m  LP  eJo 


\r  x  B ( \r) dX 


+  e<j,(r)  +  e  r  •  ET(xr)dx 

J  o 


(4.1) 


I  *• 


where  <p(r)  is  the  instantaneous  coulomb  potential  due  to 
the  nucleus  and  any  external  charges.  This  Hamiltonian 
depends  on  the  electric  and  magnetic  fields  instead  of  the 
potentials.  Expanding  the  integrals  and  the  potential,  4>(r)  , 

results  in  the  convential  multipole  expansion.  Note  that  the 
expansion  of  the  potential  contains  the  longitudinal  electric 
field,  El  ;  thus,  the  Hamiltonian  depends  on  the  total 
electric  field. 

It  is  the  intent  here  to  show  justification  for  selec¬ 
ting  Eq  (4.1)  as  the  energy  operator  based  on  the  classical 
energy  arguments  and  the  requirements  presented  in  the  pre¬ 
vious  chapter.  Detailed  calculations  for  the  full  expansion 
are  not  included,  since  this  thesis  deals  solely  with  the 
electric  dipole  portion  for  simplicity.  The  method  can  be 
extended  to  additional  multipole  terms  as  necessary. 
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Even  though  two  Hamiltonians  are  unitarily  related,  it  is 


conceivable  that  just  one  represents  the  energy  operator. 

The  interpretation  of  dynamical  variables  changes  under  a 
transformation  (19) ,  so  the  two  Hamiltonians  can  mean  very 
different  things.  Recall  that  under  a  unitary  transformation, 
the  Hamiltonian  transforms  as 


r 

* 


A  AAA  _i 

H  "  =  S  H  S 


ih  S 


(4.2) 


On  the  other  hand,  the  energy  operator  transforms  as 


A  A  A  A  i 

E'  =  S  E  S 


(4.3) 


a  %• 

a 


i. 


If  a  Hamiltonian  is  the  energy  operator  in  one  representation, 
the  transformation  of  that  Hamiltonian,  under  Eq  (4.2),  can 
represent  something  much  different  in  another  representation. 
This  allows  selection  of  Eq  (4.1)  as  the  energy  operator 
despite  the  fact  that  it  is  unitarily  related  to  Eq  (2.8), 
which  has  been  shown  not  to  accurately  represent  the  true 
energy  of  the  atom. 

For  the  derivation  of  the  electric  dipole  form,  consider 
the  transformation  of  Eq  (2.8)  with  the  electric  dipole 
approximation,  where  the  vector  potential  A  is  a  function 
of  t  only.  The  full  Hamiltonian  is 


e 

c 


A  ( t ) 


+  eO(r) 
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--  -r  -■>  ' 


H 


(4.4) 


In  Eq  (4.2),  let  the  unitary  transformation  be 


S  =  exp 


”  A 

-  i  r  • 

he 


A  ( t ) 


Using  the  fact  that  S  commutes  with  4>(r)  and  A(t)  ,  and 


that 


P,S_  = 


-  ih  li¬ 


the  electric  dipole  form  of  the  Hamiltonian  becomes 


fr  +  e<(,(r)  -  er  •  E. 


where  et  =  is  t^ie  transverse  electric  field  in  coulomb 

gauge.  Expanding  4>(r)  ,  as  in  Eqs  (3.13)  and  (3.14), 


~  *  -2  - 
H  =  2i  +  e^o  "  er  *  (EL  +  V 


where  is  the  coulomb  potential  of  the  nucleus.  This  is 

the  dipole  form  of  the  multipole  expansion;  it  contains  both 
components  of  the  electric  field.  Choosing  Eq  (4.8)  as  the 
energy  operator  is  a  different  case  than  selecting  the  full 
atomic  Hamiltonian,  Eq  (4.4). 

To  justify  use  of  the  multipole  form  of  the  Hamiltonian, 
consider  the  classical  definition  of  energy,  the  sum  of 
kinetic  and  potential  energies.  The  kinetic  energy  can  be 


written  as  pV2m  ;  in  this  case,  p  =  mv  .  The  potential 
energy  of  the  electron  depends  on  the  coulomb  potential  of 
the  nucleus,  and  the  electromagnetic  field: 


V(r)  =  e4>0  +  VE  +  VB 


(4.9) 


The  potential  energy  of  an  electron  in  an  electric  field  is 
given  by  the  expansion  (20:101) 


_  _  1  V'  V'  ^Ej 

VE  =  e<Mo)  -  er  -E(o)  -  ?L  ^Qij1¥~'~  + 


(4.10) 


These  are  the  monopole,  dipole,  and  quadrapole  terms  of  the 
electric  multipole  expansion.  is  the  quadrapole  moment 

tensor.  Similarly,  the  magnetic  multipole  energy  can  be  ex¬ 
panded  (20:148-150).  For  the  magnetic  dipole,  in  , 


VQ  =  -  m  •  B  +  ... 


(4.11) 


Thus,  the  total  energy  can  be  written  as 


E  =  + 


e  A  -  er  •  E  - 


-  m  .  B  + 


(4.12) 


This  is  of  the  same  form  of  the  multipole  Hamiltonian,  when 
operators  are  used  in  place  of  the  physical  quantities. 


Based  on  these  arguments,  the  Hamiltonian  of  Eg  (4.8)  is 


a  reasonable  choice  for  the  energy  operator  in  the  electric 
dipole  approximation.  The  electric  dipole  approximation  is 
valid  when  the  electric  field  has  a  negligible  variation  over 
the  extent  of  the  atom.  If  electric  quadrapole  or  magnetic 
dipole  interactions  are  important,  then  these  terms  should 
also  be  included  in  the  Hamiltonian. 

A  further  justification  exists  in  the  classical  limit  of 
the  equation  of  motion,  that  results  from  this  Hamiltonian. 


A 


(4.13a) 


A  > 

3H  = 

3*0 

A 

e  — s— 

3Xi 

ax. 

e  (r-E) 


In  the  dipole  approximation,  E  =  E(t)  ,  so 


(4,13b) 


(4.14) 


which  is  the  equation  of  motion  for  an  electric  field.  Note 
that  the  gradient  of  the  nuclear  potential  merely  adds  an 
additional  electric  field.  If  the  full  multipole  expansion 
were  used  as  the  Hamiltonian,  the  proper  Lorentz  force  law 


would  result. 


Gauge  Invariance  and  Stark  Effect 

The  energy  operator,  given  by  Eq  (4.8),  satisfies  both 
the  requirements  of  gauge  invariance  and  a  Stark  shift  for 
the  total  electric  field.  The  Hamiltonian  consists  totally 
of  gauge  invariant  quantities.  The  momentum  operator  repre¬ 
sents  the  kinetic  momentum.  The  potential  is  the  coulomb 
potential  of  the  nucleus.  Finally,  the  electric  dipole  inter¬ 
action  term  depends  on  the  electric  field,  which  is  gauge 
invariant.  In  another  gauge  the  Hamiltonian  becomes 


A  —2  _ 

=  2S  +  e*o  '  er 


E  =  H 


u 


(4.15) 


since  p'  =  p  and  r'  =  r  .  This  differs  from  the  previous 

A  /V  A 

case  when  p  transformed  to  p  -  eyG  ,  and  p  was  the 

0  — 

canonical  momentum  given  by  p  =  mv  +  —  A 

In  the  interaction  term,  the  electric  field  is  the  total 
field;  thus.  Stark  shifts  occur  equally  for  either  a  longitu¬ 
dinal  or  transverse  field.  This  prevents  the  violation  of 
special  relativity,  as  in  the  previous  case.  Furthermore,  the 
Stark  shift  for  a  transverse  field  provides  a  possible  means 
of  experiemental  confirmation  of  this  Hamiltonian,  if  strong 
enough  fields  can  be  generated.  The  previous  choices  pre¬ 
dicted  no  Stark  shift  for  a  transverse  electric  field,  only 
transitions  between  eigenstates.  The  other  terms  of  the 
multipole  expansion  allow  shifts  due  to  magnetic  fields  (the 
Zeeman  effect),  quadrapole  fields,  and  so  on. 


With  the  new  Hamiltonian  expression  for  the  energy  opera¬ 
tor,  quantum  descriptions  of  physical  processes  differ  from 
the  predictions  of  the  previous  two  cases.  One  significant 
difference  is  that  the  energy  operator  is  now  a  function  of 
time,  and  can  not  be  transformed  to  a  time  independent  form. 
The  eigenvalues  will  then  be  functions  of  time  also.  In  the 
previous  cases,  the  energy  operator  was  either  the  unperturbed 
Hamiltonian,  Hq  ,  or  transformed  to  the  form  of  Hq  ,  so 
the  energy  eigenvalues  were  constant.  In  the  end,  the  dif¬ 
ferences  in  the  predictions  of  the  three  Hamiltonians  provide 
a  means  to  determine  the  correct  energy  operator  by  correla¬ 
tion  with  experimental  observations. 

Dynamics  of  the  Multipole  Hamiltonian 

For  a  time  dependent  Hamiltonia,  the  solution  of  the 
Schrodinger  equation  is  not  a  stationary  state.  The  expansion 
coefficients  of  the  state  vector  are  time  dependent.  When  the 
multipole  Hamiltonian  is  considered  as  the  atomic  energy  oper¬ 
ator,  the  energy  eigenstates  and  eigenvalues  are  also  func¬ 
tions  of  time.  The  resulting  solution  of  the  Schrodinger 
equation  involves  a  state  vector  whose  expansion  coefficients 
and  basis  states  are  time  dependent.  In  other  words,  the 
interaction  with  light  causes  both  a  shift  in  the  energy 
states  and  transitions  between  the  states. 

In  the  electric  dipole  case,  the  Hamiltonian  is 

A 

/\  -“2  A  A 

H  ( t)  =  +  e<ji0(r)  -  er  •  E(t)  (4.16) 
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For  optical  frequencies  and  lower,  the  electric  dipole  form  is 
sufficient  since  there  is  little  variation  of  the  electric 


field  over  the  extent  of  the  atom.  The  magnitude  of  the  mag¬ 
netic  dipole  term  is  a  factor  of  1/c  smaller,  so  it  can  be 

A 

ignored.  Because  H(t)  is  the  energy  operator,  the  energy 
eigenstates  and  eigenvalues  satisfy  the  instantaneous  time 
independent  Schrodinger  equation  (6:289), 


H(t) | n ( t) >  =  Eft(t) | n ( t) > 


(4.17) 


The  state  vector  is  then  expressed  in  terms  of  the  instan¬ 
taneous  solutions 


4>>  =  ]T]a  (t)  |  n  ( t)  >  exp  [  En(t')dt' 
n  L  Jo 


(4.18) 


where  aR  are  time  dependent  coefficients.  The  state  vector 
satisfies  the  time  dependent  Schrodinger  equation 


ih  -2^  =  H  ( t)  |  4>> 


(4.19) 


For  ease  of  notation,  the  explicit  time  dependence  of  the 
various  terms  will  be  dropped.  Substituting  Eq  (4.18)  into 
(4.19)  and  simplifying, 


(4.20) 


Multiplying  both  sides  by  <k| 


(En-Ek)dt 


(4.21) 


This  is  an  expression  for  the  time  rate  of  change  of  the 
expansion  coefficients.  An  expression  for  <k|n>  comes  from 
differentiating  Eq  (4.17)  and  multiplying  by  <k | 


<kif? 


n>  +  <kjH|n>  =  <k 


3E 

Tt 


n 


n>  +  <k|En|n> 


(4.22) 


For  k  n  , 


3E  sE 

<k|j^  |  n>  =  <k|n>  =  0  (4.23) 

<k|H|n> =  <k |Ek|n>  (4.24) 


Thus,  Eq  (4.22)  becomes 

<k  l"at  !n> 

<  k | n>  =  ■  jj -  n  /  k  (4.25) 

Ln  ^k 

For  the  case  of  k  =  n  ,  it  is  possible  to  put  a  condition 
on  the  phase  of  the  eigenstates,  such  that  <n|n>  =  0  (6:290) 

The  final  relation  for  the  coefficients  is 


<ki!f  I 


n> 


(4.26) 


i 


where 


“nk(t|  = 


En(t)-Elc(t) 


(4.27) 


This  relation  describes  transitions  between  the  time  depen¬ 
dent  energy  eigenstates. 

The  form  of  Eq  (4.26)  is  considerably  different  than  its 
counterparts  for  the  first  two  Hamiltonians,  Eqs  (2.19)  and 
(2.31).  The  transition  matrix  depends  on  3H/9t  ,  which  is 

-  er*DE/3t  in  the  electric  dipole  form.  The  exponential 
term  contains  an  integral,  since  the  transition  frequency  is 
time  dependent.  Finally,  recall  that  the  eigenstates  are 
functions  of  time,  as  well.  In  the  previous  case  of  the 
full  Hamiltonian,  Eq  (2.30),  the  eigenstates  are  functions  of 
time,  but  are  unitarily  related  to  the  unperturbed  eigen¬ 
states,  simplifying  the  solution.  Unfortunately,  this  is  not 
the  case  for  Eq  (4.26).  These  differences  translate  to  a 
significantly  more  difficult  solution  when  applying  Eq  (4.26) 
to  a  particular  experimental  setup. 

Physical  predictions  using  the  multipole  form  for  the 
energy  operator  should  differ  from  those  of  the  unperturbed 
or  full  atomic  Hamiltonian.  The  choice  of  the  energy  opera¬ 
tor  is  most  important  in  describing  such  phenomenon  as  multi¬ 
photon  processes,  atomic  lineshapes,  non-exponential  decay 
rates,  and  absorption  sidebands.  The  difference  in  predic¬ 
tions  of  the  three  models  is  a  means  of  confirming  the 
correct  choice  for  the  energy  operator.  Already,  a  possi- 
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bility  exists  in  the  Stark  effect  for  a  transverse  electric 
field.  The  multipole  Hamiltonian  is  the  only  one  that  has  an 
equal  dependence  on  the  longitudinal  and  transverse  electric 
fields.  Measuring  a  Stark  shift  for  a  transverse  field 
should  distinguish  between  the  various  Hamiltonians.  This 
is  not  a  trivial  task,  since  very  strong  fields  are  required. 
Another  possibility  is  in  the  explanation  of  the  absorption 
sidebands  for  OCS  in  a  time  varying  electric  field,  observed 
by  Townes  et .al . (21 : 273-279 ;  22).  Townes  uses  a  perturbation 
calculation,  but  the  multipole  Hamiltonian  may  yield  a  better 
explanation.  This  is  an  area  for  further  study. 

In  any  case,  the  multipole  Hamiltonian  is  a  plausible 
choice  for  the  energy  operator,  and  could  be  used  to  explain 
many  physical  phenomena.  In  the  next  chapter,  a  multiphoton 
process  is  examined  to  show  the  different  predictions  of  the 
three  energy  operators,  and  provide  a  means  of  experimental 


conf irmation. 


V. 


Three  Photon  Absorption  Calculations 

The  ultimate  confirmation  of  the  multipole  form  of  the 
Hamiltonian  as  the  energy  operator  depends  on  correlation 
with  experimental  results.  Comparison  between  the  predictions 
of  Eq  (4.26)  and  those  of  the  other  two  choices,  Eqs  (2.20) 
and  (2.31),  should  indicate  the  proper  choice.  One  physical 
process  which  has  a  significant  difference  between  predic¬ 
tions  is  three  photon  absorption.  This  last  chapter  details 
the  calculations  for  three  photon  absorption  of  the  three 
Hamiltonians.  Order  of  magnitude  differences  exist,  as  well 
as  a  difference  in  the  resonant  frequencies. 

To  simplify  the  calculations,  a  two  level  system  is 
assumed.  This  is  a  reasonable  assumption,  since  many  atomic 
systems  can  be  approximated  in  this  way  or  forced  into  such 
a  condition.  In  addition,  spontaneous  emission  is  ignored 
in  the  calculations,  which  allows  much  simplification  of  the 
rate  equations.  If  the  experimental  evaluation  is  made  on  a 
time  scale  short  compared  to  the  spontaneous  emission  rate, 
the  assumption  is  valid. 

/\  A 

Calculations  for  H  and 

The  Schrodinger  equation  for  the  full  atomic  Hamiltonian, 
H  ,  can  be  expressed  in  the  form  of  Eq  (2.31), 


.  1  . 
a  =  -rz-  >  <  e  |  -er »E  e_>a_e 
n  in  i—J  n 1  'mm 


iuj  t 
nm 


m 


(5.1) 
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where  the  coefficients  an  represent  the  expansion  coeffici¬ 
ents  of  the  true  energy  eigenstates.  The  eigenstates  in  the 
transition  matrix  are  those  of  the  unperturbed  Hamiltonian, 


£n> 


(5.2) 


If  the  atom  is  assumed  initially  to  be  in  the  ground  state  of 
a  two  level  system,  then  the  time  evolution  of  a2  describes 
the  absorption  rate  for  a  given  electric  field.  The  solution 
of  Eq  (5.1)  involves  the  use  of  time  dependent  perturbation 
theory.  For  three  photon  absorption,  third  order  perturba¬ 
tion  terms  are  required. 

For  a  two  level  system,  Eq  (5.1)  becomes  a  pair  of 
coupled  rate  equations: 


-iwofc 


al  "  ih<£l  -er*E|e2>a2e 


(5.3a) 


1  —  —  +iuvt 

a2  =  Ih<£2i'er*El£l>ale 


(5.3b) 


where  the  transition  frequency  is  defined  as 


OJ 


o 


(5.4) 


and  H '  =  -er*E  is  the  perturbation.  Let  the  electric  field 
be  linearly  polarized,  and  defined  as  E  =  ^  cosut  ?\  ,  where 


is  the  orientation  unit  vector.  Also  define 


(5.5) 


U  =  e  <  £]_  [r  I  ^2  >  *  ^ 


The  transition  matrix  for  this  two  level  system  is 


H  '  = 

nm 


o  -^ycosoJt 


-^ycosoot  o 


Eqs  (5.3)  then  become 


a.^  =  if^coscat  e 


-id)  t 
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a2  =  ifta^cosdit  e 


id)  t 
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where  o  =  ^ 


ft  =  ,  is  the  Rabi  frequency. 

From  time  dependent  perturbation  theory,  let 


H  '  =  -  Xer  «E 


an(t) 


t)  =  ai0)  (t)  +  Xa  (1)  (t)  +  \2i 


(2) 


(t)  +  Xa'  '  (t)  +  X  an  (t)  + 
n  n  n 


(5.6) 


(5.7a) 


(5.7b) 


(5.8) 


(5.9) 


Substituting  into  Eq  (5.1)  and  equating  powers  of  x  , 


(m) 
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m-l)  ,  I*,  ^nk^  dt  ' 
<k  H  n>e 


(5.10) 


This  equation  expresses  the  mth  order  perturbation  in  terms 
of  the  previous  one.  Eqs  (5.7)  then  become 
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= 
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(5.11b) 


If  the  two  level  system  is  initially  in  the  ground  state, 

(t)  =  a  (0)  =  6  i  .  The  first  order  perturbation  is 
n  n  nl 

aj1^  =  0  ,  and 
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(5.12) 


This  is  the  familiar  first  order  perturbation  calculation, 
which  describes  resonant  single  photon  absorption.  When  to 
approaches  the  transition  frequency,  tu0  ,  the  second  term 
becomes  relatively  large,  and  takes  the  form  of  a  sine 

function.  The  sine  function  has  a  series  of  sidelobes  which 
are  unphysical.  This  results  from  the  unnatural  assumption 
of  an  instantaneous  switching  on  of  the  electric  field,  and 
is  not  physically  significant. 

The  second  order  perturbation  is 
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The  third  order  perturbation  uses  Eq  (5.14) 
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Expanding  the  integrand, 
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This  equation  finally  yields  the  three  photon  absorption 


term.  When  cj=toQ/3  ,  the  last  term  becomes  slowly  varying, 
and  thus  increases  greatly  after  integration.  It  is  the  only 
physically  significant  term,  although  other  terms  of  the  form 
I/oj^+oj  are  also  slowly  varying.  These  additional  terms  are 
essentially  constant  over  the  range  of  interest,  a)Q/3  ; 

thus,  it  is  felt  they  are  unphysical  and  a  result  of  a  limita¬ 
tion  in  the  perturbation  approach.  Further  justif ication  for 
ignoring  these  terms  is  given  when  the  final  solution  is  com¬ 
pared  to  an  exact  solution  of  this  problem.  Dropping  the 
unnecessary  terms. 
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1  1  2 


i  (ojQ-3to)  t 
e _ 
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dt 


i  _ /lyiit 
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(5.17) 


For  three  photon  absorption,  u-uq/3  ,  there  are  no 
contributing  factors  from  the  first  order  perturbation, 

(3 ) 

Eq  (5.12).  The  upper  level  coefficient  is  then  a2  =  a2 
The  probability  of  being  in  the  upper  energy  state  is 


P2  (t)  =  |  a2  ( t )  |  2 


(5.18) 


which  can  be  written  as 


(5.19) 


p  (t)  =  •&— 

256  u  4 


sin  j(a)0-3(ij)  t 
(u>0”3u))  2 


This  equation  represents  the  probability  for  three  photon 
absorption  for  the  full  atomic  Hamiltonian.  The  resonant 
frequency  is  u  =  cjq/3 

As  a  comparison,  the  exact  solution  from  Shirley  (23:983- 
985)  is 


P2(t) 
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(5.20) 


The  two  equations  are  very  close  for  small  t.  When 
2  4  4 

(3u)-u)o)  >>  2  /4ajQ  ,  the  equations  are  the  same,  and  as  u> 
approaches  uq/3u  ,  they  have  the  same  limit.  If  t  is 
large,  the  perturbation  method  is  not  valid,  since  the  per¬ 
turbation  is  assumed  to  be  small.  The  favorable  comparison 
gives  further  justification  for  ignoring  the  extra  terms  of 
Eq  (5.16).  The  results  from  the  perturbation  calculation  are 
also  comparable  to  the  method  of  Autler  and  Townes  (22) . 

Treating  the  unperturbed  Hamiltonian  as  the  energy 
operator  results  in  a  perturbation  of  the  form 

/v 

ti'  =  -e/inc  A’p  ,  and  after  Eq  (2.19),  the  Schrodinger  equa¬ 
tion  becomes 


-  Sr  A-pi 
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iU)  t 
nm 


(5.21) 


The  only  difference  between  Eq  (5.21)  and  Eq  (5.1)  is  the 
transition  matrix.  Kobe  (13:626)  shows  that  the  matrix 
elements  are  related  such  that 


<£n! me  A*P| 


Em>  =  “|f  <£n  jer,E  |  £m> 


(5.2 


Thus,  the  development  tor  the  unperturbed  Hamiltonian  paral¬ 
lels  the  previous  calculation  except  for  a  constant.  Each 
perturbation  term  changes  by  a  factor  of  uQ/ui  .  For  a 

third  order  perturbation,  the  expansion  coefficient  is  multi 

3 

plied  by  ( co0/ to )  •  The  probability  of  being  in  the  upper 

energy  state  is  now 


P 


(t) 


(5.2 


Thus,  for  Hq  ,  the  probability  for  three  photon  absorption 
is 


(t)  =  9" 


'A«p 


256  4 
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sin  ^ (w0-3cj)  t 

(coQ-3co)2 


(5.2 


The  result  is  about  three  orders  of  magntiude  greater  than 
that  of  the  full  Hamiltonian. 


Multipole  Hamiltonian  Calculation 

The  final  calculation  is  for  the  multipole  Hamiltonian, 
which  presents  some  distinct  differences,  primarily  due  to 
the  time  dependence  of  the  energy  eigenstates  and  eigenvalues 


The  general  method  of  solution  is  to  solve  for  the  instan¬ 
taneous  eigenstates  and  eigenvalues  of  Eq  (4.17),  treating  t 
as  a  constant  parameter.  The  instantaneous  solutions  can  then 
be  used  in  Eq  (4.26),  which  can  then  be  solved  using  time 
dependent  perturbation  theory. 

For  a  two  level  system,  Eq  (4.26)  becomes  a  pair  of 
coupled  rate  equations 

a  -  _ifta)21dt" 

='5^7  <M  if  i2> e  °  (5-25a) 

al  afi  ^o'"21’^ 

J2  -  HJjr  <2lff  I  L>  e  <5-25b> 

Solution  of  these  equations  is  not  trivial  due  to  the  time 
dependence  of  the  eigenstates  and  eigenvalues.  The  energy 
eigenstates,  j n>  ,  are  the  instantaneous  solutions  of  the 
time  independent  Schrodinger  equation 

H(t) | n ( t) >  =  En(t) | n (t) >  (5.26) 

The  solution  follows  the  method  of  Cohen-Tannoud j i  (1:406-408) 
420-423)  .  The  multipole  Hamiltonian  can  be  written  as 


Using  the  unperturbed  eigenstates  |  en>  as  the  representation 

A 

basis  set,  H  can  be  expressed  in  matrix  form: 


H  =  <  e  H  e  >  = 
nm  n 1  1  m 


COSCOt 


-U^coscot 


(5.28) 


where  is  defined  by  Eq  (5.5).  Diagonalizing  Hnm  deter¬ 

mines  the  true  energy  eigenvalues. 
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(5.29a) 

(5.29b) 


where  again 
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The  instantaneous  transition  frequency  is 
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(5.30) 


which  oscillates  in  time,  above  the  unperturbed  value.  The 
oscillation  may  be  the  primary  cause  for  such  phenomenon  as 
multiphoton  absorption  and  absorption  sidebands.  The  eigen¬ 
states,  in  the  representation  are 
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i>  =  I  2  >  = 

\-sine/  \cos0  / 


(5.31) 


where 
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(5.32) 


The  transition  matrix  of  Eq  (5.25)  can  now  be  determined 
using  Eqs  (5.5)  and  (5.31): 
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From  Eq  (5.11) , 
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Thus,  Eqs  (5.25)  become 
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(5.36b) 


These  are  exact  equations,  but  approximations  must  be  used  to 
arrive  at  manageable  solutions. 


From  time  dependent  perturbation  theory,  let 
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a_ (t)  =  a(0)(t)  +  Xa(1)(t)  +  X2  a<2) (t)  +  ...  (5.38) 

n  n  n  n 


Substituting  into  Eq  (4.26)  and  equating  equal  powers  of  X 
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This  equation  is  comparable  to  Eq  (5.10). 

If  the  two  level  system  is  initially  in  the  ground  state, 
the  first  order  perturbation  for  the  upper  energy  level  is 


a21} (t) 
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(5.40) 


The  instantaneous  transition  frequency  co '  can  be  expanded 
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(5.41) 
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Higher  order  terms  can  be  ignored  provided  Q  <<  1 


The  integral  in  the  exponential  becomes 
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>  2 

where  =  ajQ  +  ft  /4u>o  is  the  shifted  transition  frequency, 
as  will  be  shown  later.  This  shift  is  a  function  of  the  field 
strength,  £  ,  which  demonstrates  the  Stark  effect  for  a 

transverse  field.  The  exponential  can  be  expanded  as  a  series 
of  Bessel  functions  (24:361): 
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For  z  <<  1  (24:358) , 
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Using  Jn  =  (-l)nJn 
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Ignoring  the  higher  order  terms,  the  exponential  becomes 
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Expanding  the  demoninator  of  Eq  (5.40), 
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Finally,  Eq  (5.40)  becomes 
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which  is  directly  integrable,  although  messy.  Expanding 

3 

integrand  and  dropping  terms  of  higher  order  than  H 
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Notice  that,  only  odd  harmonics  of  u>  are  present  in  the  expo¬ 
nential  terms.  Absorption  occurs  only  when  a  term  becomes 
slowly  varying,  when  the  exponential  argument  approaches  zero. 
Therefore,  only  an  odd  number  of  photons  can  be  absorbed,  in 
this  model.  The  conclusion  corresponds  to  the  conservation 
of  angular  momentum.  The  angular  momentum  must  change  by  one 
unit  in  the  transition,  and  photons  have  one  unit  of  angular 
momentum.  Only  odd  combinations  are  possible. 

Integrating  Eq  (5.49) 
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gets  relatively  large  whenever  w  approaches  , 

the  shifted  transition  frequency,  or  uJ/3  ,  three  photon 
absorption.  If  01=03^/3  ,  then  only  the  last  term  is  of  im¬ 
portance,  and  since  ^q~^q  • 
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It  is  necessary  to  calculate  higher  order  perturbation 

3 

terms  since  they  contribute  terms  of  order  .  Recall 

that  since  the  atom  is  initially  in  the  ground  state,  the 
perturbation  terms  ai^  =  ai^  = 
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But 
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which  is  similar  in  form  to  Eq  (5.40).  Expanding  the  inte¬ 
grand  as  before,  using  Eq  (5.50)  for  a2^  •  anc^  dropping 

higher  order  terms  of  ft  results  in  an  integrable  form.  The 
solution  of  a^  is  then  used  in  the  calculation  of  the 

third  order  perturbation: 
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Again,  if  u)=:u>J/3  ,  only  one  physically  significant  term  in 

the  integrand  is  slowly  varying.  Other  unphysical  terms  are 
also  slowly  varying,  but  these  can  be  ignored  as  in  the 
previous  case,  Eq  (5.16).  The  perturbation  term  becomes 
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The  result  of  the  integration  can  be  expressed  as 
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The  final  expression  for  a2  is  a2  =  a2^  +  a2^ 
from  Eqs  (5.51)  and  (5.55)  is 
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The  probability  of 

being  in  the  upper  energy  level  is 
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which  represents  the  transition  for  three  photon  absorption 
with  the  multipole  Hamiltonian.  The  resonant  frequency  for 
this  transition  is 


ST 

12  uj 

♦  o 


The  results  for  the  multipole  Hamiltonian,  Eqs  (5.57) 
and  (5.58)  are  much  different  than  those  of  the  other  two 
Hamiltonians . 


Comparison  of  Predictions 


Calculations  for  the  three  Hamiltonians  shows  consider¬ 
able  differences  in  the  predictions  for  the  probability  of  a 
two  level  atom,  initially  in  the  ground  state,  to  transition 
to  the  upper  energy  level.  In  each  case,  time  dependent 
perturbation  theory  was  used  to  determine  the  probability. 
The  method  is  valid  for  relatively  short  time  periods, 
otherwise  effects  such  as  Rabi  oscillations  and  spontaneous 


emission  become  important.  In  addition  to  the  transition 
rate  difference,  the  multipole  Hamiltonian  case  also  showed 
a  Stark  shift  for  a  transverse  electric  field.  The  first  two 


Hamiltonians  do  not  account  for  such  an  effect. 

The  differences  in  the  upper  level  predictions  is  most 
easily  seen  for  radiation  that  is  exactly  on  resonance, 
u;  =  (jjq/3  .  The  probabilities  for  the  three  cases  are 
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Thus , 


P2Ap  ~  10  P2rE  ~  2x10  P2MH 


(5.62) 


This  demonstrates  the  order  of  magnitude  differences  between 
the  three  choices.  These  differences  are  directly  translated 
to  comparable  differences  in  the  absolute  transition  rates  or 
the  Einstein  B  coefficients.  The  rate  for  the  multipole 
Hamiltonian  is  the  smallest  by  a  factor  of  twenty.  If  the 
absolute  rate  could  be  measured,  it  would  strongly  indicate 
the  correct  choice  for  the  energy  operator. 


Absolute  transition  rates  for  a  three  photon  process 
are,  unfortunately  for  this  problem,  difficult  to  measure. 
First,  there  is  a  requirement  for  a  very  strong  field  to 
obtain  an  appreciable  amount  of  transitions.  Assuming  a 
sufficient  source  is  available,  there  are  problems  with  the 
media  itself.  In  crystals,  the  problem  of  wave  vector 
matching  (25:76-78)  severely  hinders  absolute  measurements. 

The  problem  is  non-exponential  absorption  due  to  inter¬ 
ference  effects  in  the  crystal.  Another  possibility  for  an 
experiment  is  to  measure  the  deflection  of  an  atomic  beam  by 
radiation  pressure.  An  incident  light  wave  can  cause  a  beam 
deflection  based  on  the  absorption  rate  (26)  .  The  primary 
difficulty  with  this  configuration  is  in  generating  a  stable 
beam  that  is  slow  enough  to  have  a  reasonable  interaction 
time . 

The  other  major  difference  in  the  predictions  was  the 
resonant  frequency  shift  predicted  only  in  the  multipole 
Hamiltonian  case.  Considering  either  the  unperturbed  or  the 
full  atomic  Hamiltonian  as  the  energy  operator  allows  no 
shift  due  to  a  transverse  electric  field,  although  such  shifts 
have  been  observed  for  time  varying  fields  (22)  .  The  predic¬ 
tion  of  this  Stark  shift  for  the  multipole  Hamiltonian  may  be 
the  best  argument  for  this  choice  as  the  energy  operator. 
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VI .  Conclusion 


The  choice  of  the  correct  form  of  the  Hamiltonian  as  the 
energy  operator  has  significant  implications  for  physical  pre¬ 
dictions  in  quantum  mechanics.  Historically,  two  choices  have 
been  predominant:  the  unperturbed  Hamiltonian  and  the  full 
atomic  Hamiltonian.  Each  results  in  different  predictions  as 
indicated  by  their  respective  transition  matrices.  This 
difference  is  often  clouded  by  the  equivalence  of  these  two 
forms  under  a  unitary  transformation.  The  argument  for 
equivalence  is  valid,  based  on  interpretation  of  the  eigen¬ 
states,  but  it  does  not  address  the  identification  of  the 
atomic  energy  operator.  This  thesis  has  argued  that  both  of 
these  choices  are  limited  due  to  conflicts  with  either  observed 
phenomenon  or  physical  laws,  and  neither  represents  the  true 
energy  of  the  atom. 

Another  possibility  for  the  energy  operator  exists  in  the 
multipole  expansion  form  of  the  Hamiltonian.  It  is  derived 
from  the  classical  multipole  energy  expansion,  and  has  the 
correct  classical  limit  of  the  Lorentz  force  law.  The  multi¬ 
pole  form  is  gauge  invariant,  and  contains  the  total  field, 
so  it  does  not  violate  special  relativity.  The  dynamics  of 
this  form  are  much  more  complicated  than  the  previous  cases, 
which  compounds  the  difficulty  of  making  physical  predictions. 
Thus,  for  the  purposes  of  first  order  approximations,  the 
unperturbed  or  full  atomic  Hamiltonian  should  prove  to  be  an 


adequate,  though  not  completely  accurate,  choice  for  the 
energy  operator.  For  finer  applications,  the  multipole  form 
may  provide  the  better  model  for  the  atom,  and  hence,  a  more 
accurate  explanation  for  physical  processes. 

Calculations  for  three  photon  absorption  demonstrate 
the  differences  between  predictions  of  the  three  choices  for 
the  energy  operator.  Order  of  magnitude  differences  exist 
for  the  probability  of  transitioning  to  the  upper  state  of 
a  two  level  atom.  These  directly  translate  to  comparable 
differences  in  the  absolute  transition  rates.  The  multipole 
Hamiltonian  has  a  further  distinction  of  predicting  a  Stark 
shift  for  a  transverse  electric  field.  The  other  two  Hamil¬ 
tonians  do  not  allow  this  shift  in  the  respective  models. 

Experimental  confirmation  of  the  absorption  prediction 
or  the  Stark  effect  due  to  a  transverse  field  is  the  final 
proof  of  this  thesis.  Unfortunately,  measurement  of  an 
absolute  three  photon  absorption  rate  is  a  difficult  task, 
and  has  not  been  performed.  This  is  an  area  for  further 
investigation.  In  addition,  Stark  shifts  have  been  observed 
in  numerous  experiments  (22) .  Further  study  of  those  experi¬ 
mental  results  and  the  methods  of  analysis  is  warranted. 

The  stark  effect  may  provide  the  best  confirmation  of  the 
correct  energy  operator. 

Another  suggestion  for  future  work  is  to  consider  the 
affect  of  spontaneous  emission  in  the  absorption  calculations 
The  steady  state  solution  to  a  density  matrix  perturbation 


calculation  should  give  the  absorption  rate  for  long  inter¬ 
action  times,  and  spurious  zero  frequency  terms,  previously 
disregarded,  should  be  damped  out. 

Other  recommendations  include  the  investigation  of 
absorption  sidebands,  observed  by  Townes  (21:273-279)',  and 
atomic  lineshapes.  These  areas  may  provide  further  arguments 
for  choosing  the  multipole  Hamiltonian  as  the  atomic  energy 


operator 
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